Testing integrability with a single bit of quantum information 
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We show that deterministic quantum computing with a single bit (DQCl) can determine whether 
the classical limit of a quantum system is chaotic or integrable using 0{N) physical resources, where 
A'^ is the dimension of the Hilbert space of the system under study. This is a square root improvement 
over all known classical procedures. Our study relies strictly on the random matrix conjecture. We 
also present numerical results for the nonlinear kicked top. 
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I. OVERVIEW 

After an initial triumph at solving mathematical prob- 
lems (see Q for an overview), a large fraction of the re- 
searches in the field of quantum information processing 
has shifted to its original motivation: the simulation of 
quantum systems fj| . It is now well established H, ^ |^ 
that the evolution produced by certain classes of Hamil- 
tonians can be simulated efficiently on a universal quan- 
tum processor. However, extracting useful information 
from the physical simulation is a problem who's complex- 
ity has been underestimated. Indeed, the ability to sim- 
ulate the dynamics of a system does not grant one with 
the ability to evaluate efficiently all physical quantities 
of interest. These quantities (e.g. spectral properties) 
are usually measured experimentally on a (exponentially) 
large number of physical systems — macroscopic sam- 
ples. A direct quantum simulation, on the other hand, 
can only reproduce the statistical output of a single quan- 
tum system which yields drastically less information than 
what is learned from costly classical simulations. Thus, 
it is not clear at this point whether quantum simulators 
can always outperform their classical analogues. 

Some of these spectral properties play a central role 
in the study of quantized chaotic systems. One partic- 
ular question of interest is whether the classical limit of 
a quantum system exhibits regular or chaotic motion. It 
has became widely accepted (see 0, Q and references 
therein) that the answer to this question is hidden in 
some spectral properties of the system, which can be re- 
produced by those of canonical random matrices with 
the appropriate symmetries. Given a description of the 
Hamiltonian of the system, the best known algorithms 
evaluating these "signatures of chaos" require classical 
computing resources which grow at least as fast as N'^, 
the square of the dimension of the Hilbert space of the 
system under study. Indeed, a close inspection of these 
algorithms show that they require either matrix multipli- 
cation, diagonalization, or evaluation of a determinant 



[g. Since such a growth is intractable on any conven- 
tional computer (remember that N grows exponentially 
with the size of the physical system) , it is quite natural 
to try to tackle this problem with a quantum computer. 
In recent years, this interest has lead to the demonstra- 
tion that the standard model of quantum computation 
can simulate efficiently the dynamics of a few quantized 
chaotic models ,9., 10] ; unfortunately, none of these 
proposals indicate how to circumvent the measurement 
problem mentioned above. 

Recent work by Emerson et al. [TJ proposes to study 
statistical properties of the system's eigenvectors rela- 
tive to a perturbation as a signature of chaos. They also 
provide an efficient procedure to measure these statis- 
tics using the standard model of quantum computation. 
Their motivation for this work was to test the validity of 
the signature. Indeed, it is not clearly established that 
this signature is universal, perturbation independent and, 
most importantly that the decay time does not scale with 
the size of the system. 

Here, we concentrate on a different model (presumably 
weaker) : deterministic quantum computation with a sin- 
gle pscudo pure bit (DQCl) which was introduced in 
In this setting, the initial state of the K + \ qubits com- 
puter is p = {i^ll-He|0)(0|} ® 27^11 where < e < 1 is 
a constant. Note that, from a computational complexity 
point of view, this is equivalent to a model where the 
state of the first qubit is pure while the other ones are 
completely random; we shall therefore assume that e = 1 
in the remaining of the paper. The final answer is given 
by a finite accuracy evaluation of the average value of ct^ 
on the first qubit. As for the dynamics, we assume that 
we are gifted with the ability to exert coherent control 
over one and two qubits at a time. This model is of par- 
ticular interest since it is weaker than the computational 
model offered by liquid state nuclear magnetic resonance 
(NMR) quantum computing 13]. Such a computing de- 
vice, we shall show, can test for integrability using 0{N) 
physical resources, given that the dynamics of the sys- 
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tern of interest is efficiently simulatable on the standard 
model of quantum computation without ancillary pure 
qubits (or, more precisely, with no more than 0(logK) 
ancillary pure qubits). 

In order to do so, we must first relate the theory under- 
lying the spectral property at the heart of our study; this 
is done in Sec.|n| We then show how it can be evaluated 
with 0{N) physical resources in the DQCl model. In 
Sec. IIVI we present numerical results for canonical ran- 
dom matrices as well as for a physical map, the nonlinear 
kicked top. Finally, we conclude with a summary of our 
results and discussion of possible extensions. 



II. LEVEL DISTRIBUTION 

In the theory of quantum chaos, a key role is played 
by the statistics of eigenvalues 0, 0- In the case of 
systems with a periodically time varying Hamiltonian 
the central dynamical object is the Floquet operator, 
F — T'[exp{— i H{t)dt}], that maps the state from one 
time to a time exactly one modulation period T later, T 
is the time ordering operator. The eigenvalues of F lie 
on the unit circle and may be parameterized in terms of 



eigenphases, or quasi-energies, as F\(j)j 
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The random matrix conjecture asserts that the statis- 
tics of eigenvalues of chaotic systems (dynamical systems 
and maps) is typically well modeled by the statistics of 
the eigenvalues of random matrices (hermitian Hamilto- 
nians and unitary Floquet operators) with the appropri- 
ate symmetries ja, Q| • While many important mathemat- 
ical results underpin the conjecture, a rigorous proof is 
lacking and support rests on a very large accumulation 
of numerical results. 

An integrable system, by definition, possesses as many 
symmetries — constant of motion — as degrees of free- 
dom. One can thus write the system's Hamiltonian as 
the direct sum of independent Hamiltonians acting on 
smaller subspaces; one for each values of the constants 
of motion. Some spectral properties of these Hamiltoni- 
ans can thus be reproduced by those of matrices that are 
the direct sum of independent random Hermitian opera- 
tors. The distribution characterizing the entire spectrum 
is therefore given by the superposition of many indepen- 
dent spectra; as a consequence, the correlations between 
levels vanish. Thus, one might expect that the near- 
est neighbors level spacing distribution (LSD) follow a 
Poisson law prob{(f>j+i — (pj — S) = P{S) ^ e~^^ , a 
straightforward consequence of their statistical indepen- 
dence. This is indeed observed experimentally, numeri- 
cally, and most importantly can be derived formally ■ 

On the other hand, chaotic systems possess no or just 
a few symmetries. It can be shown that the LSD — 
aside from the systematic degeneracy following the sym- 
metries — obeys a power law ^(5*) ~ S^e'"^ . The 
parameter /3 characterizes the symmetries of the system; 
it is equal to 1 when the system possesses a time reversal 
symmetry and some geometric invariance, 2 when it has 



no symmetries, and 4 when it has a time reversal sym- 
metry with Kramer's degeneracy. Similarly, we will refer 
to the Poisson ensemble — the characteristic ensemble of 
integrable systems — as /3 = 0. 

The exact form of the LSD is not relevant to us; 
we shall capitalize on the crucial distinct behavior of 
P{S — !■ 0) for chaotic and regular systems. In the former 
case, P{S) reaches a minimum at 5* = 0: the levels tend 
to repel each other. In the latter case, P{S) is maximal 
at S* = 0, a consequence of the levels statistical indepen- 
dence called clustering. 

With these considerations, one can predict the behav- 
ior of the ensemble average form factors 
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from which most spectral properties can be extracted. 
For regular systems, Tr{F} = J2j e"*"^^ behaves like 
the end point of a random walk in the complex plane: 
each step having unit length and uncorrelated random 
orientation ipj. After N steps, the average distance 
from the origin is expected to be so we should find 
Ti = N. For times n > 1, the analysis is identical; if 
the angles {(f>j} are statistically independent, then so are 

{(/)^"^ = n(/)jmod(27r)}, n taking positive integer values. 
We conclude that the ensemble average form factors of 
integrable systems should be time independent and equal 
to N. 

For chaotic systems, more elaborate calculations are 
required for the ensemble average form factors. They can 
be found in ; here we shall simply give an approximate 
result for < n < iV (accuracy of order 10"^) known as 
the Wigncr surmises: 
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Although simple arguments could not have indicated 
the exact behavior of these form factors, we could have 
guessed their general form: they are initially very small 
Ti <^ N, and, as n grows, they reach the same value as 
the Poisson ensemble. Here, Tr{F} is analog to a anti- 
correlated random walk in the complex plane composed 
of N unit steps. As a consequence of level repulsion, 
each steps tend to be oriented in different directions; the 
probability of finding two steps oriented within an angle 
e decreases as e^^^ . Thus, the distance from the origin 
after N of these anti-correlated steps should definitely be 
smaller than \/N which is the expected value for uncor- 
related steps. As n grows, the phases mod (27r) get 
wrapped around the unit circle; the effect is analogue to 
superposing n independent spectral distributions, blur- 
ring out the correlations. When n ^ N, one should thus 
expect a behavior similar to the Poisson ensemble. 

It should be noted that the few symmetries of a chaotic 
system may slightly affect the predictions of Eq. |21 The 
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average r„ were evaluated for fixed values of the con- 
stant of motion. In what follows, we shall often neglect 
this point for sake of simplicity. Nevertheless as long 
as the number of invariant subspaces is small (<C VN) 
this omission will not affect our conclusions. For ex- 
ample, if a chaotic system possesses a symmetry which 
breaks its Hilbert space into k equal invariant subspaces, 
the small n behavior T„ n will be transformed into 
T„ ~ k^n <C N, which is all that really matters to us. 
One can circumvent this issue when some exact symme- 
tries of the system are known: it suffices to simulate the 
dynamics of the system within an invariant subspace. 

In the light of this analysis, it may seem that form 
factors constitute a powerful tool to distinguish between 
classically regular and chaotic systems. In particular, r„ 
should clearly identify each regime for small values of 
n. Nevertheless, the form factor r„ of a fixed Floquct 
operator F will generally fluctuate about the ensemble 
average T„. Thus, we seek a signature of an ensemble 
property on a single element drawn from this ensemble. 

The solution is to use a version of the ergodic theorem. 
If we normalize out the explicit time dependence of the 
form factors, an average over a time interval An repro- 
duces the effect of an ensemble average. More precisely, 
one can show 01 that 
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converges to 1 with a variance bounded by 1/An. For 
large N , we can thus use the first An <C N form factors 
to determine whether the Floquet operator belongs to a 
polynomial or Poisson ensemble. Since the value of r„ 
— hence the matrix ensemble — are needed to compute 
Eq. 13 we shall proceed by hypothesis testing: for which 
choice of r„ (T„ = N regular, T„ ~ n chaotic) does Eq.El 
converge to 1? In other words, we need to determine 
which of the two variables 
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is most probably drawn from a distribution centered at 
1 with l/VAn standard deviation. If we restrict our at- 
tention to a regime where An <C iV, both hypothesis 
cannot have high probabilities simultaneously On 
the other hand, when the probabilities of both hypothesis 
are low, the test is inconclusive. Nevertheless, remember 
that the presence of symmetries in a chaotic system shifts 
the value of the distribution by a factor k^ where k is the 
number of invariant subspaces. For fc^ ^ N , this should 
be clearly distinguishable from the value of a regular sys- 
tem. This should not be seen as a bug but a feature of 
our approach allowing one to estimate k, the number of 
invariant subspaces. 

Applying this test to a particular dynamical system 
would require one to compute the spectrum of the Flo- 



quet operator. If one were to try and simulate a dynam- 
ical map on a quantum computer with K qubits, a di- 
rect computation would require determining all = 2^ 
eigenvalues. In the next section we will construct a quan- 
tum circuit that would enable the form factors themselves 
to be extracted with 0{N) physical resources thus allow- 
ing a direct test of non-integrability that circumvented 
the need to explicitly compute all eigenvalues. 



III. QUANTUM ALGORITHM 

The DQCl algorithm evaluating the form factor is 
based on the idea reported in 15] of using a quantum 
computer as a spectrometer. The circuit is shown at 
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FIG. 1: Quantum circuit evaluating the trace of . The 
gates B!g are rotation in the Bloch sphere by an angle 6 around 
axis k — X or y. When k is set to x, we get the real part of 
the trace while k = y yields the imaginary part. 

to efficiently simulate the dynamics of the system under 
study so the gate i^" only requires a polynomial (in n 
and K) number of elementary gates to be constructed. 
Here, it is not i^" we wish to implement but a coherently 
controlled version of it, i.e. a linear gate acting on AT 4- 1 
qubits which applies F" to the last K qubits when the 
first qubit is in state |1) and doesn't do anything when 
it is in state |0). Given the circuit for F", standard tech- 
niques can be used to construct a controlled version of it 
at polynomial cost [T^ . 

It should also be emphasized that the K qubits on 
which the Floquet operator is applied generate a Hilbert 
space of dimension 2^ which might be larger than the 
simulated system's Hilbert space. Thus, when applying 
F to those qubits, one really applies F (B U where ide- 
ally U is the identity operator on 2^ — A^ states; it can 
be any other unitary operator as long as its trace can 
be evaluated. The effect of these extra dimensions will 
be to add a contribution Tr{U}/N to the output signal 
which should be systematically subtracted as we shall 
henceforth assume. 

The output of this computation will be the real and 
imaginary part of {Tr{F''^})/2^ when the last rotation 
is made about axis k = x and k = y respectively. Thus, 
our task is to distinguish between a signal whose ampli- 
tude is of order 1 /N (chaotic dynamics) and one of order 
\/^/N (regular dynamics) which can be achieved using 
0{N) physical resources. In the special case of NMR 
quantum computing, one can for example increase the 
size of the sample by a factor A^ as the size of the system 
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increases, or simply repeat the procedure N times and 
sum up the outputs. We thus get a quadratic advantage 
over all known classical algorithms. 



IV. NUMERICAL RESULTS 

A. Random Matrices 

Before applying our general proposal to a physical 
model, we give a numerical example illustrating the main 
results used from random matrix theory: the ergodic the- 
orem of Eq. ^ In order to estimate the average and vari- 
ance of to and ti in a given universal matrix ensemble, 
we draw many random matrices J/''^'^ from the ensemble 
and numerically evaluate each quantity. As an example, 
we have generated 50 random matrices from the (3 = 2 
ensemble — the set of unitary matrices with no symme- 
tries. This is illustrated on Fig.|21 where the matrices are 
of size 600 x 600. For each random matrix C/'^'^) drawn 
from this ensemble, we can compute ti{U^''^) as functions 
of An. Two such curves (dashed) are plotted on Fig. |21 
By applying this procedure to many samples (here 50), 
we can estimate the average of ti and its fluctuations: 

-.50 .50 



The average (ti) and mean deviation a ~ 
V ^ S'i'e also plotted on Fig. |3 (heavy 

and light full line respectively): as expected, ti con- 
verges to 1 as 1/V An. The same procedure can be 
applied to to; nevertheless, since ti does converge to 1 
in this ensemble, to obviously does not since it differs 
by a factor of roughly N/An ~ 20 for the range of 
An we have studied. Of course, this difference would 
vanish when An approaches N since the form factor of 
any universal ensemble converge to those of the Poisson 
ensemble (see Sec. m; this is why we must restrict our 
study to An ^ N. The same conclusions can be reached 
for the other ensembles characterizing chaotic systems, 
i.e. /3 =1, 2, and 4. 

Similarly, had the matrices J/'-'"'^ been drawn from the 
f3 = ensemble — the set of matrices characterizing reg- 
ular systems — we would have observed to converging 
to 1 as l/VAn while ti, smaller by a factor of roughly 
An/iV, would roughly vanished. From these considera- 
tions, the hypothesis test "to converges to 1" versus "ti 
converges to 1" allows us to discriminate between random 
matrices drawn from /? = and those drawn from one of 
the P = 1, 2, or 3, with a probability of error decreasing 
as 1/VAn. Thus, as long as the random matrix conjec- 
ture holds, it should also allow to discriminate between 
regular and chaotic motion. 
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FIG. 2: The two dash lines show ti as a function of An (Eq.|lll 
for two random unitary matrices drawn from the ensemble 
/3 = 2. The heavy full line is the value of ti averaged over 50 
such random matrices while the light line shows its variance, 
which drops as 1/ V An as expected. 

B. Kicked Top 

We now focus our attention on a physical model of 
great interest for its good agreement with random matrix 
theory: the nonlinear kicked top. We write the Floquet 
operator in its most general form following Haake 0], 
F = U^UyU^ with 

Uk = cxp I ~^ 2j'+l ~ ^'^'"^'^ I 

where the Jk, k = x, y, and z, are the canonical angular 
momentum operators. We conveniently define a param- 
eter vector p = [ux, ay,az,Tx,Ty,Tz). Some authors use 
a restricted form of this Floquet operator where only 
and ay are non-zero. Since [F, J^] = 0, the value of the 
angular momentum j — which appears in Eq. — is 
conserved. The dimension of the Hilbert space is simply 
given by = 2j + 1. 

By adequately choosing the parameters p, the kicked 
top can be either in a regular or chaotic regime, see 
for more details. Thus, we can evaluate to and ti of Eq.^ 
in both regimes and verify that they indeed allow to dis- 
criminate between them. This is presented on Figs. Inland 
2] for different values of the total angular momentum j. 
On Fig. |31 the system is in a regular regime; we have only 
plotted to since ti is larger by a factor proportional to j so 
clearly does not converge to 1. Similarly, only the value 
of ti is exhibited of Fig. 0] Notice that while the ergodic 
averaging decreases the fluctuations, it is not essential to 
discriminate between regular and chaotic. Indeed, the 
scale of the fluctuation is extremely small compared to j 
which is the factor by which to and ti differ. 

The analogy with a random walk in the plane can also 
be illustrated graphically. On Fig. Owe have plotted the 
sum of the eigenvalues vectorially. The apparent struc- 
ture of the vectors is purely artificial, the eigenphases 
were ordered in increasing order (the sum of vectors is 
obviously a commutative operation); we have chosen this 
ordering to facilitate the presentation. 
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FIG. 3: Value of to (Eq. of the kicked top in a regular 
regime pr — (0,0,1,0,0,10) as in for different values of 
j. Dashed curve: An = 1 so it is simply Tr{F} /N . Full 
curve: To decrease the fluctuation, we have used the ergodic 
averaging over the first An = 30 normalized form factors 
Tr{F"}/N, n = 1,2, ...30. 




j 

FIG. 4: Value of ii (Eq. 0J of the kicked top in a chaotic 
regime pc = (1.1, 1, 1, 4, 0, 10) as in Q for different values of j. 
Dashed curve: An = 1 so it is simply Tr{F}. Full curve: To 
decrease the fluctuation, we have used the ergodic averaging 
over the first An — 30 normalized form factors Tr{F"} /n, 
n = 1,2, ...30. 



The effect of LSD are striking on Fig. [S| The light 
vectors (chaotic regime) are arranged in an almost perfect 
circle; eigenphases tend to be equally separated. On the 
other hand, the heavy vectors (regular regime) are quite 
often aligned in an almost straight line; a manifestation of 
level clustering. As a consequence, the heavy vectors end 
up further apart from the origin than do the light vectors; 
on average, these distances differ by a factor ^/N. 

Finally, we can use the form factor to study the transi- 
tion between regular and chaotic motion. To do so, we let 
the parameter vector continuously vary from its regular 
value to its chaotic value: p = (1 — e)pr + epc (see figure 
captions O and 0J. For e = 0, the expected value of to 
is 1. As e increases, the system enters a chaotic regime; 
when chaos has fully developed, should vanish as 1/iV. 
This is indeed observed on Fig. |^ where we have plotted 
^0 as a function of e for different system sizes. Moreover, 
the results indicate that the transition to chaos becomes 
more sensible as the size of the system increases. 
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FIG. 6: Value of to (Eq.|4|l with An = 30 for different system 
size: Full heavy line j — 50; Dashed line j — 100; Light full 
line j = 200. 
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FIG. 5: Vectorial representation of eigenphases: 
J {cos sin where the 0^ have been ordered in 
increasing order. The a and r parameters of the Floquet 
operator Eq. Elare tuned so the system is in a regular regime 
(heavy vectors) and a chaotic regime (light vectors) as in 
Figs. |3] and 2] The value of j is 20 so each curve contains 41 
vectors. 



V. CONCLUSION 

We have shown that, using a single bit of quantum 
information, we can test whether the spectrum of a uni- 
tary matrix obeys a Poisson or polynomial law. Under 
the random matrix conjecture, this can be used to de- 
termine whether the system has a regular or chaotic be- 
havior in its classical limit. The idea relies on estimat- 
ing the averaged form factor using the ergodic theorem 
which roughly states that a time average can reproduce 
an ensemble average. The form factors in a regular and 
chaotic regime differ by a factor of ^/N and the output 
signal of our computation decreases as 1/N: the required 
physical resources thus scale as TV. This is a quadratic 
improvement over all known classical algorithms. We are 
currently investigating a different signature of quantum 
chaos which might not suffer from this signal loss and 
hence, could offer an exponential speed up. 

This result gives a new insight on the nature of the po- 
tential computational speed up offered by quantum me- 
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chanics. In particular, it provides a strong argument to- 
wards the computational power of mixed states quantum 
computing. 
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